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ABSTRACT 

We present four new seasons of optical monitoring data and six epochs of X-ray photometry for the 
doubly-imaged lensed quasar Q J0158-4325. The high-amplitude, short-period microlensing variability 
for which this system is known has historically precluded a time delay measurement by conventional 
methods. We attempt to circumvent this limitation by application of a Monte Carlo microlensing 
analysis technique, but we are only able to prove that the delay must have the expected sign (im- 
age A leads image B). Despite our failure to robustly measure the time delay, we successfully model 
the microlensing at optical and X-ray wavelengths to find a half light radius for soft X-ray emis- 
sion \og(r 1 / 2 .x,soft/ cm ) = 14-31q5, an u PP er limit on the half-light radius for hard X-ray emission 
log(r 1 / 2 .x,/iarrf/ cm ) < 14-6 and a refined estimate of the inclination-corrected scale radius of the 
optical i?-band (rest frame 3100 A) continuum emission region of log(r s /cm) = 15.6 ± 0.3. 

Subject headings: cosmology: observations — accretion, accretion disks — dark matter — gravitational 
lensing — quasars: general 



1. INTRODUCTION 

Quasars can be highly variable. Since this variability 
arises in the quasar continuum source itself, it appears 
in all the images of a gravitationally lensed quasar, but 
at different times in each image due to differences in 
the gravitational potential of the lens galaxy and the 
geometric length of the paths of each image (Refsdal 
19641: iCooke fc Kantows ki 1975). As first pointed out by 
Chang fc R efsdal (1979;), the stars near each lensed im- 
age also cause temporal variations in the magnification of 
each image, leading to additional variability. This "mi- 
crolensing" variability is not correlated between the im- 
ages, so it functions as both a difficult source of noise in 
time delay measurements and as a probe of quasar struc- 
ture because the size of the emission region ultimately 
controls the amplitude of the microlensing variability. 

Traditional time delay measurements require either 
that the microlensing variability is of sufficiently small 
amplitude so as to be insignificant in the analy- 
sis or that its timescale is significantly longer than 
that of the quasar's intrinsic variability, allowing for 
some sort of parametric modeling and subtraction (e.g. 
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Ofek fc Mao3[200l IParaficz et al.ll2006l IKochanek et al.l 
20061 : IPoindexter et al.l 120081 ). Attempting to measure 
time delays in systems with significant short-timescale 
microlensing is both difficult and error prone (see 
lEigenbrod et al.ll2005l for quantitative characterizations 
of the influence o f micro lensing on time delay estimates). 
In iMorgan et al.l (|2008| ) , we applied a new approach to 
analyzing the light curves of two lensed quasars with 
significant microlensing and time variability, in which 
we attempted to simultaneously model the microlens- 
ing vari ability usin g the Bayesian Monte Carlo tech- 
nique of [Kochanek (2004). In our new approach, we 
marginalized over the physically realistic models of the 
microlensing signal to yield a measurement of the most 
likely time delay. We successfully confirmed the time 
delay of HE 1104-1805 but failed to measure the delay 
in Q J0158-4325. In the present paper, we supplement 
our previously published light curves with 4 new seasons 
of monitoring data from two observatories and make a 
second attempt. 

In the models of Morgan et al.l ()2008[ ). we could also 
marginalize over the time delay uncertainties and use the 
microlensing to study the structure of the quasar. The 
amplitudes and the timescales of the variability due to 
quasar microlensing are a function of a number of physi- 
cal variables such as the continuum source size, the mean 
properties of the lens in the vicinity of the lensed images, 
and the relative velocities of the source, the lens and the 
observer. During the last decade a variety of techniques 
have been developed to analyze the effects of microlens- 
ing on lensed quasar images. Single-epoch photometry in 
which one or more of the quasar image fluxes are anoma- 
lous (i.e. deviates significantly from the predictions of the 
macroscopic lens model or shows significant wavelength 
dependencies) can provide (prior dependent) estimates 
of source sizes and const rain temper ature profiles (e.g. 
Poolev et all [2001 [20071 iBate et all [200l iFlovd et all 
20091: iBlackburne et all 120111 iMosquera et all \20T\i 
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Figure 1. Top Panel: QJ0158-4325 it-band light curves for the magnitudes of images A and B + 0.1. Data from the SMARTS 1.3m 
are plotted as black squares. Data from the 1.2m Euler Swiss Telscope are plotted as red triangles. Observations failing a seeing limit of 
FWHM < 2"0 are not shown. This eliminated 20 points from the light curves. 



iMediavilla et al.ll2011l ). These studies must also assume 
a mean mass for the microlensing stars. 

Stronger and less prior-dependent constraints can be 
obtained from analyses of quasar microlensin g variabil- 
ity. We have employed the technique of IKochanekl 
( 2004) to analyze uncorrelated variability in quasar opti- 
cal lightcurves to yield a correlation between the sizes of 
quasar accretion disk s and the masses of their black holes 
(|Morgan et al.l [2010h . meas urements of the structu re of 
quasar accretion disks (e.g. IPoindexter et al.ll2008T) . and 
even a constraint on the orien tation of a quasar's ac- 
cretio n disk (in Q2237+0305, IPoindexter fc Koc hanck 
120101) . We have also modified the IKochanekl (|2004D 
technique to simultaneously analyze optical and X-ray 
lightcurves, yielding measurements of the size of the 
X-ray emitting regions in four systems, PG1115-I-080, 
HE 1104-1805, RXJ 1131-1231 and HE 0435-1223 
([Morgan et alJl200llChartas et alj|2009t iDai et alllMol: 
iBlackburne et al.l 120121 respectively). In those papers, 
we found that the X-ray continuum emission emerges 
from region ~ 10f g in radius, where r g is the black 
hole's gravitational radius r g = GMbh / c 2 ■ This is 
at least an order of magnitude smaller than sizes mea- 
su red in the UV /optic al, consistent with the arguments 
of iPoolev et al.l (|2007l ). These results seem to disfavor 
X-ray conti nuum emission from an extended disk-co rona 
model (e.g. IHaardt fc Maraschl 19911: IMerlonill2003|) , but 
more measurements are needed both to confirm these 
results and to begin searching for correlations between 
the X-ray emitting regions and other physical proper- 
ties of the quasars. In this paper, we add to the sample 
of quasars studied at X-ray wavelengths by applying our 
Monte Carlo microlensing analysis technique to the com- 
bined optical and X-ray lightcurves of QJ0158-4325 to 
yield size estimates for the hard and soft X-ray contin- 
uum emission regions. Making use of our new constraints 
on the time delay, we are also able to significantly refine 
our estimate of the optical accretion disk size. 

In §[2] we begin by descr ibing our Chandra X-ray obser- 
vations ([Chen et al.ll2012l) . followed by a presentation of 
our new optical monitoring data and a brief review of our 
data reduction process. In §[3J we describe our time delay 
measurement attempts and in § 0] we present the results 
from our Monte Carlo microlensing simulations. In § [5l 



Table 1 

Q J0158-4325 Light Curves 



HJD 


Seeing 


QSO A (mags) 


QSO B (mags) 


Source 


4316.305 


1.55 


1.869 ± 0.012 


2.541 ±0.018 


Euler 


4329.389 


1.34 


1.881 ± 0.008 


2.550 ±0.013 


SMARTS 


4329.396 


1.61 


1.842 ± 0.009 


2.558 ±0.015 


Euler 


4335.299 


2.10 


(1.872 ±0.009) 


(2.567 ±0.015) 


SMARTS 


4342.326 


1.31 


1.833 ± 0.014 


2.589 ±0.025 


Euler 


4347.368 


1.80 


1.861 ±0.010 


2.531 ±0.016 


Euler 


4348.295 


1.20 


1.853 ± 0.008 


2.565 ±0.012 


SMARTS 


4354.319 


2.01 


(1.876 ± 0.010) 


(2.544 ±0.015) 


Euler 


4357.293 


1.30 


1.861 ± 0.008 


2.593 ±0.012 


SMARTS 



Note. — HJD is the Heliocentric Julian Day - 2450000 days. The 
QSO A&B columns give the magnitudes of the quasar images relative 
to the comparison stars. A few points in the light curves with seeing 
FWHM > 2'.'0 (in parentheses) were not used in the analysis. Table 1 
will be published in its entirety in the electronic edition of The Astro- 
physical Journal. A portion is shown here for guidance regarding its 
form and content. 

we conclude with a discussion of the implications for the 
physics of the optical and X-ray continuum sources in 
this quasar. All calculations in this paper assume a flat 
ACDM cosmology with h = 0.7, tt M = 0.3 and fl\ = 0.7. 

2. OBSERVATIONS AND LIGHT CURVE REDUCTION 

We imaged the d oubly lensed quasar Q J0158-4325 
(jMorgan et al.lll999T) in X-rays usin g the ACIS imaging 



spectrometer (jGarmire et aL I2003D on the Chan dra X- 



Ray Observatory ( Chandra, IWcisskop f et al.ll2002D on six 
occasions between 2009 November 6 and 2010 October 6. 
These observations were a component of a larger Chan- 
dra Cycle 11 m onitoring program , the details of which 
are published in iChen et all (|2012t ). We divided the full 
observed frame X-ray spectral range (0.4-8.0 keV) into 
soft (0.4-1.3 keV) and hard (1.3-8.0 keV) energy bands. 
We analyze the soft and hard X-ray variability separately 
in order to constrain the energy structure of the X-ray 
continuum source. 

We have monitored Q J0158-4325 for eight seasons in 
the i?-band using the SMARTS 1.3 m telescope with th e 
ANDICAM optical/infrared camera (|DePov et al.l2003f FI 
and using the 1.2m Euler Swiss Telescope as a part of 

8 http:/ /www. astronomy.ohio-state.edu/ANDICAM/ 
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Figure 2. Example of Monte Carlo fits to the observed flux ratios in the lensed quasar Q J0158— 4325. Top Panel: Six epochs of the 
observed A/B X-ray flux ratio from our Chandra monitoring program shown as mx A — rax B - An example of a good Monte Carlo fit to 
these measurements is also shown. Bottom Panel: The Q J0158-4325 optical difference light curve for a trial time delay of At^B = — 8 days. 
Data are shown with error bars. The curve is a good fit (\ 2 /N^of ~ 1-2) from our Monte Carlo analysis. Note its failure to fully match all 
features of the trial difference light curve. 



the COSMOGRAIL3 project. We analyzed the optical 
images from both observatories si multaneously, employ - 
ing the deconvolution algorithm of iMagain et al.l £l998). 
In this technique, we measure the point spread function 
(PSF) in each frame using the same 5 reference stars in 
the field. The astrometry of the lens components is held 
fixed to the values measured in Hubble Space Telescope 
iJ-band images from the CfA- Arizona Space Telescope 
Lens Sur vey (CASTLE S ), the details of which were pub- 
lished in [Morgan ct al. (2008). Using the corresponding 
PSFs, we simultaneously fit the photometry of the point 
sources and a regularized image of the lensing galaxy to 
all the frames. The resulting QSO image fluxes are cal- 
ibrated relative to the deconvolved field stars. Finally, 
we apply a small multiplicative and additive correction 
to minimize the dispersion in the overlapping sections of 

9 http://www.cosmograil.org/ 



the Euler and SMARTS lightcurves. This empirical ad- 
justment is required to compensate for small color terms 
introduced by differences in the filters and quantum ef- 
ficiency curves of the CCDs. The new photometric data 
are presented in Table [TJ and our composite eight-season 
light curves are displayed in Figure [TJ 

3. ANALYSIS 

3.1. Lightcurve Preparation 

We began our analysis by shifting the image A light 
curve by a reasonable set of trial time delays (—30 days 
< AtAB — tA ~ ts < +30 days) separated by one day 
intervals, where A leads B for AtAB < 0. The model 
prediction for a singular isothermal sphere (SIS) lens 
model is that A leads B by AtAB — —12.4 days. We 
chose to shift image A because it is less variable than 
image B. Our analysis algorithm requires contemporane- 
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Figure 3. Left: Probability distribution for the time delay AtAB = tA ~ tfl m Q J0158-4325, where we restricted the extrapolation of 
the light curves to be less than 5 days. Negative delays are now strongly favored over positive delays as predicted by lens models, but our 
analysis did not converge on a single value for the delay. The solid vertical line at AtAB = —12.4 days represents the predicted delay of 
a SIS lens model. Right: Results of Legendre polynomial fits to the observed lightcurves; relative probability as a function of time delay 
AtAB = tA — is is plotted. The intrinsic and microlensing variability are fit with 30 th and 3 rd order polynomials, respectively. Note the 
absence of a single, well-defined peak. 



ous pairs of A/B data points, in most cases necessitat- 
ing linear interpolation between bracketing epochs. For 
shifted points that fell into the inter-season gaps, we per- 
mitted 5 days of linear extrapolation, and we increased 
the size of the measurement uncertainties on the extrap- 
olated points using our estimate of the quasar's intrinsic 
variability structure function. Since the Bayesian inte- 
gral shown in Eq. Q] is actually a discrete sum in which 
the statistical weight of each trial time delay must be 
held fixed, all trial lightcurves must have the same num- 
ber of epochs, forcing us to discard some data points 
for the shorter trial delays . The details of this selection 
technique are described in lMorgan et al.1 (|2008l ). 

We also generated a second set of trial light curves, 
permitting unlimited extrapolation in the inter-season 
gaps using t h e Damped Random Wa lk ( DRW) models of 
lKelh^(|20ll ; lKozlowski eTaLl (|2010h an dlMacLeod et al " 



Using the implementation of iKozlowsk i et al 



2010| ). we first modeled the raw lightcurves, yielding 



a long term variability structure function measurement 
SFqo = 0.337 mags and timescale r = 1867 days for 
image A and SFqo = 0.280 mags and timescale r = 
419.3 days for image B. These models also provide the 
structure function estimate used in our first approach. 
The obvious advantage of this treatment is that no data 
must be discarded, however the DRW process models the 
intrinsic variability in the light curves of only one image 
at a time such that the variability in the DRW points is 
not correlated with the contemporary points in the other 
image. In the absence of the strong microlensing signal, 
we would m odel the light cur ves jointly, following the 
procedure of iPress et al.l (|1992P ). The introduction of ad- 
ditional uncorrelated variability can be significant, par- 
ticularly for the longer trial delays, forcing the models 
to fit the differences in the DRW interpolations as mi- 
crolensing. This further blurred any effects of the time 
delay, so we do not report on these results. 



Because we have only 6 epochs of X-ray observations, 
we are unable to shift the X-ray measurements by the 
trial time delay, since doing so would eliminate most of 
our data points. Both our macroscopic mass models and 
this present work predict a relatively short time delay 
— 30 < AtAB < days, but it is likely that some in- 
trinsic variability also contaminates the X-ray flux ra- 
tios. To account for this, we estimated the variability 
structure function of the (sparsely sampled) X-ray light 
curves and then conservatively assumed a 30 day time 
delay, yielding an additional systematic uncertainty of 
0.1 magnitudes in the X-ray flux ratios. This additional 
systematic uncertainty limited the precision of our Monte 
Carlo analysis, but it cannot be avoided until the time 
delay is measured. 

3.2. Simultaneous Microlensing and Time Delay 
Analysis 

Using the same techn ique that we described in detail 
in I Morgan et al.l ((2008), we appl ied the M o nte C arlo 
microlensing analysis method of iKochanekl (|2004t ) to 
the light curves from each trial time delay, in essence 
attempting to reproduce the observed lightcurves by 
random combinations of the unknown variables: the 
lens galaxy stellar mass fraction, the median microlens 
mass (M) and the effective velocity between source, 
lens and observer v e . W e used the recentl y measured 
lens redshift z t = 0.317 (|Faure et al.l I2009D to update 
our set of ten lens models, each consisti ng of concen- 
tric NFW JN avarro. Frenk fc Wh4ti[l996l ) and de Vau- 
couleurs mass profiles. We parametrize the stellar mass 
fraction 0.1 < Jm/l < 1 in our model sequence by com- 
paring the fractional mass of the de Vaucouleurs compo- 
nent to the mass of a constant M/L model (Jm/l = 1); 
varying Jm/l between 0.1 and 1.0 in uniform steps. For 
each model, we generate four unique sets of microlens- 
ing magnification patterns using the P 3 M method de- 
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scribed in Kochanck ( 20041). The magnification patterns 
are 4096 x 4096 pixel arrays with an outer scale radius of 
20R E , where R E = 3.4 x 10 16 ((A/)/0.3Af Q ) 1 / 2 is the 
Einstein radius, yielding a source plane pixel scale of 
1.7 x 10 14 ((M)/0.3M Q ) 1 / 2 cm. We performed four sep- 
arate Monte Carlo simulations, one per set of magnifi- 
cation patterns. In each simulation, we attempted 10 s 
trials per mass model, and in each trial we attempt to fit 
all 61 trial time delays, for a grand total of 2.4 x 10 11 tri- 
als. We discard all solutions in which x 2 /^dof > 5.0, 
since they contribute negligibly to the Bayesian integrals 
and consume unnecessary disk space. An example of a 
set of good fits to the observed flux ratios is shown in 
Figured 

We use Bayes' theorem to marginalize over the results 
from our Monte Carlo simulation such that the probabil- 
ity of a given time delay At 



P(At\D)<x J P(D\p, At)P(p)P(At)dp, 



(1) 



where P(D\p, At) is the probability of fitting the obser- 
vations in a given trial, P(p) sets the priors on the mi- 
crolensing variables (see lKochane k 2004; Kochanc k et al.l 
120071) and P(At) is the (uniform) prior on the time delay. 

The probability density for the time delay is displayed 
in Figure [3] To first order, our analysis appears to have 
failed, since it did not converge on a robust, single- valued 
estimate for the time delay. We did succeed, however, 
in effectively eliminating the possibility of positive val- 
ues for the delay At AB, which is an imp rovement over 
our earlier results in IMorgan et al.l ((2008). Nonetheless, 
the present probability density for the delay is (at least) 
double-valued. We also attempt ed a traditional polyno - 
mial fit time delay analysis (e.g. iKochanek et alj (2006). 
in which we fit the quasar's intrinsic and microlensing 
variability using Legendre polynomials. The best-fit re- 
sults from this attempt, based on a 30 th order fit to the 
intrinsic and 3 rd order fit to the microlensing variabil- 
ity, are also displayed in Figure Again, as was the 
case with the Monte Carlo technique, our polynomial fits 
did not converge on a single value for the delay, in fact 
the (inferior) polynomial fitting technique still allows for 
some significant probability of a positive (B leading A) 
time delay. Nevertheless, the Monte Carlo procedure al- 
lowed us to study the structure of the quasar with proper 
statistical weighting of the uncertainties created by the 
time delay, ultimately enabling a successful joint fit to the 
observed X-ray and optical lightcurves using our Monte 
Carlo method. 

For the simultaneous optical and X-ray anal y sis, w e 
employed the technique detailed in iDai et all (|2010f) . 
in which a large number of good fits to the opti- 
cal light curves are generated first. To enable bet- 
ter resolution at smaller angular scales, we used a set 
of 8192 x 8192 magnification patterns with an outer 
scale of IORe, providing a source plane pixel scale of 
4.2 x 10 13 ((M) /0.3M Q ) 1/2 cm, ~ 2 times larger than the 
black hole's gravitational radius at r g — 2.5 x 10 13 cm. 
We discuss the limitations of this pixel scale in § |U We 
also tested the Monte Carlo simulation using the original 
20Re magnification patterns to confirm that the IORe 
patterns from the low /m/l mass models did not intro- 
duce any significant systematic errors. We saved all phys- 



ical variables from the good optical fits (e.g. trajectory, 
effective transverse velocity v e ) and then fit the X-ray 
light curves using the same path across the magnifica- 
tion pattern but a different range of source sizes. An 
example of a good fit from this analysis is also displayed 
in Fig. [2] We then analyze the combined X-ray and op- 
tical solutions to yield a set of joint probability densities 
for the variables of interest, namely the X-ray and opti- 
cal source sizes. Since the time delay was unknown, we 
performed this analysis for a range of assumed delays, 
but the X-ray size was insensitive to the resulting minor 
changes in the optical solutions. We also tested our opti- 
cal solutions for the possibility of 20 — 30% line emission 
contamination from the quasar's un-microlensed broad 
line region, but t his led to only modest changes in th e 
size estimates (see IMorgan et alJl2010HDai et alll feOlO). 

Similar to our analysis in IMorgan et al.l (|2008f ). the un- 
known time delay introduces short-timescale noise in the 
optical difference lightcurve that the Monte Carlo code 
can successfully reproduce using non-physical transverse 
velocities. To eliminate this non-physicality in our mod- 
els, we applied a prior on the effective transverse veloc- 
ity 10 (M/M ) 1 / 2 kms- 1 < v e < 10 4 (Af/M Q ) 1 / 2 kms- 1 
where v e is in Einstein units. To convert all variables 
from Einstein units of (M /Mq) 1 / 2 into physical units, 
we combined the resulting probability density for the ef- 
fective transverse velocity P(v e ) with the distribution 
for the true transverse velocity P{v e ). We generated 
the physical mo del for P(v e ) following the method of 
IKocha nek (2004), where we estimated the peculiar ve- 
lo cities of the source and lens u sing the power-law fits 
of iMosquera fc Kochanekl (|201 If ) . We found the trans- 
verse velocity of the observer by projection onto to the 
CMB dipole vcmb — 328.7 km s _1 , and we used a pe c,i = 
277km s^ 1 and <J pec , s — 248km s^ 1 for the peculiar veloc- 
ities of lens and source, respectively. We estimated the 
velocity dispersion in the lens galaxy cr* = 203 kms _1 
using the monopole moment of a singular isothermal el- 
lipsoid (SIE) lens model. 

4. RESULTS: OPTICAL AND X-RAY CONTINUUM SIZES 
AND LENS GALAXY PROPERTIES 

In Figure IK we display the probability densities for 
the size of the continuum emission region at optical and 
X-ray wavelengths. In our microlensing simulations, we 
assume a simple thin accretion disk model where the disk 
radia tes locally as a blackbody (e.g. iShakura fc Sunvaevl 
Il973h . Since we expect the optical accretion disk to be 
many times larger than the black hole's gravitational ra- 
dius r g — GMbh /c 2 , we assume that we can safely ignore 
the inner edge of the accre tion disk and the as sociated 
relativistic corrections (see M organ et al.l l20Toh . With 
these assumptions, the surface brightness at rest wave- 
length A rest is 
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where the scale length 
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Figure 4. Joint probability distributions for the X-ray half-light radius Ti/2,x (Top Panel) and inclination-corrected optical source size 
r a (Bottom Panel) of Q J0158— 4325, assuming cosi = 0.5. In the X-ray panel, the solid, dotted, and dashed curves show the probability 
densities for the half light radius of the full band (rest frame 0.9-18.4 keV), soft band (rest frame 0.9-3.0 keV) and hard band (rest frame 
3.0-18.4 keV) continuum source, respectively. In the optical panel, the solid curve shows the estimate for r s from the simult aneous analysis 
of the X-ray and optical light curves. The dashed curve displays our original estimate of r B from Morgan ct al. (2008) based on only 
4 seasons of optical monitoring data. The vertical lines show the last stable orbit in the Schwarzschild metric at 6r g and at lr g for a 



maximally rotating black hole in the Kerr metric, 
the Mg II emission line width. 



The black hole mass M BH = 1.6 X 10 s M was estimated bv lPene et al.i 1 120061 ) using 



45G\$ est M BH M 
I6n 6 h n c 2 



1/3 



(3) 



= 9.7 x 10 1 
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V 10 9 Mr. 
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( L 



1/3 



cm 



is the radius at which the disk temperature matches the 
wavelength, kT\ rtsgt = h p c/\ res t, h p is the Planck con- 
stant, k is the Boltzmann constant, Mbh is the black 
hole mass, M is the mass accretion rate, L/Le is the 
luminosity in units of the Eddington luminosity, and 
r] = L/(Mc 2 ) is the accretion efficiency. In previous 
investigations (e.g. iMorgan et al.ll2010t iPoindexter et al.l 



2008) we found accretion disk temperature profiles con- 
sistent with the T(r) oc r~ 3 / 4 predictions of thin disk 
theory, so we report the size of the optical continuum 
emission region in terms of this model. 

Microlensing observations can estimate the tempera- 
ture slope, since for T a r"* 3 the disk size should scale 
as r\ oc A 1 /' 3 = A 4 / 3 for standard thin disk theory. Al- 
though there is general agreement that hotter regions 
(i.e. shorter wavelengths) are more compact, observa- 
tional estimates of the slope (3 have signifi cant uncertain- 
ties a n d dispersion between s t udies (e.g. lAnguita et all 



2008t lEigenbrod et all l2008t IPoindexter et all 120081; 
Bate et al l 120081: iFlovd et al.l l2009t iBlackburne et all 
20111 ). FortunatelvT lMortonson et al.l (|2005l ) determined 
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Figure 5. Top Panel: Probability density for the ratio of 
the size of the soft X-ray continuum emission region to the 
size of the hard X-ray emission region. The expectation value 
\°g{ri/2,x,soft/ri/2,x,hard) = 0-2t°;? is consistent with unity, but 
tends toward a larger soft X-ray source. The solid vertical line at a 
size ratio of unity is plotted to guide the eye. Bottom Panel: Prob- 
ability density for the ratio of the optical half-light radius r"i/2 j0 pt 
to the half light radius of the full-band X-ray source T\ji t x- The 
most likely ratio is log(r 1 /2,opt/ r l/2,Js:) = 1-7 ± 0.5. 

that microlensing statistics are insensitive to the choice 
of model surface brightness profile and depend primar- 
ily on the half-light radius. Thus, the size we report 
here can easily be compared to other estimates using the 
half-light radius ri/ 2 = 2.44r s . Also, since microlensing 
statistics are actually sensitive to the projected area of 
the quasar's accretion disk, we must account for the un- 
known inclination angle i. We assume that the optical 
emission comes from a geometrically thin disk so that the 
de-projected estimate of the size is r s = r Sy0 b s / '(cos i) 1 ' 12 
for random orientations. For simplicity, we use the same 
surface brightness distribution in X-rays but assume that 
the emission is geometrically thick, so we report the half- 
light radius r l / 2X - 

In iMorgan et all (|2010T) we demonstrated that in the 
context of the thin accretion disk model, a quasar's spe- 
cific luminosity L\ can also be used to estimate the size 
of its accretion disk at that rest-frame wavelength. In the 
observed-frame /-band, for example, the so-called "flux- 
based" thin disk scale radius is 



r s ./ = 2.83 x 10 



15 



1 



A/. 



obs 



VcosI 

.3/2 



D. 



OS 



10 -0.2(7-19) h - 



(4) 



There is now significant evidence in the literature 
that flux(luminosity)-based size estimates are con- 
sistently s maller than microlensing-based measure- 
ments (e.g. iPoindexter et al.ll200 8: Hainl ine et all 120121 ; 
IBlackburne et all 12012ft . and in IMorgan et all (|2010f ). 
we found this to be true in 10 out of the 11 quasars 
in our sample. The one outlier in that sample was 
Q J0158-4325, whose flux- and microlensing-based sizes 
were consistent. Integration of our updated r s proba- 



bility density and correction for an average inclination 
of cosi = 0.5 yields an optical accretion disk size es- 
timate of log[(r s /cm)(cos(i)/0.5) 1 / 2 ] = 15.6 ± 0.3 at 
Xrest = 0.31/im, the rest-frame center of the i?-band. So 
in our new, more precise measurement, the flux-based 
size (\og(r s ji ux /cm) = 15.2 ± 0.1) is now ~ 0.4 dex 
smaller than the microlensing-based size, consistent with 
the discrepancies in the other systems. 

Probability densities for the half-light radii of the full, 
hard and soft X-ray continuum emission regions are also 
displayed in Figure 2] Since there is no simple spatial 
model for the X-ray continuum source, we report the X- 
ray sizes in terms of the half-light radius ?~i/2,x for X-ray 
emission at rest-frame energies 0.9-18.4 keV (full band), 
0.9 - 3.0 keV (soft band) and 3.0-18.4 keV (hard band). 
Integration of the probability density for the combined 
hard and soft X-ray flux yields an estimate for the full 
band half light radius log(r 1/2)X /cm) = 14.2±° j. The 
soft and hard size estimates are consistent with that of 
the full band and each other, where the half light radius 
for soft X-ray emission is log(ri/2,x, so/t/ cm ) = 14.3±q 5 
and the upper limit on the half-light radius for hard X- 
ray emission is log(riy 2 .x, /mrd/ cm ) < 14.6. Because the 
pixel scale in our magnification patterns is 7.6 x 10 13 cm 
for a 1 Mq star, the lower limit on our X-ray size es- 
timates is not fully resolved. In essence, once a trial 
source becomes much smaller than the pixel scale, all so- 
lutions become equally likely. In the case of the hard 
X-ray source, this prevented any convergence on a lower 
limit, so we also compiled a probability density for the 
ratio between the soft and hard X-ray half-light radii 



log(ri 



1/2, X,soft/ r l/2,X, hard) 



0.2 



+0.6 
0.7- 



This probability 



density is plotted in Figure [5] along with the ratio of the 
optical half-light radius to that in the full X-ray band, 

which peaks at login / 2 , opt /r 1 / 2 , x ) = 1.7 ± .5. 

In other systems. IMorgan et al.l (|2008l ): iDai et all 
(pOlofl and IBlackburne et all (|2012D have all found X- 
ray continuum regions that seem to be concentrated in 
the vicinity of the black hole's last stable orbit. In 
Q J015 8-4325, the blac k hole mass M BH = 1.6 x 
10 s Mq (jPeng et alll2006l ) implies a gravitational radius 
log(r 9 /cm) = 13.4 and last stable orbit log(r sta (,/ e /cm) = 
14.2 at 6r g in the Schwarzschild metric. Consistent 
with our results in HE 0435-1223, RXJ 1131-1231 and 
PG1115+080, we find that the majority of the X-ray 
emission in Q J0158-4325 emerges from a region in the 
vicinity of the last stable orbit. We also find that the 
soft and hard X-ray continuum emission regions seem 
to be the same size, consistent with the findings of 
IBlackburne eTaTI ((20121) for HE 0435-1223. 

In Figure we display the results of marginalizing 
the fits from our Monte Carlo simulation over all vari- 
ables other than the stellar mass fraction in the macro- 
scopic mass model. The resulting constraint on fuiL 
favors a dark-matter dominated halo, which seems to be 
marginally consistent with our expectations for a time 
delay of approximately —8 days. 

5. DISCUSSION AND CONCLUSION 

To demonstrate the improvement in accuracy gained 
by adding four seasons of new optical monitoring data 
and simultaneous analysis with X-ray lightcurves, we also 
show the size distribution from our first analysis of this 
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Figure 6. Probability distribution for /m/l> the fraction of the 
lens galaxy mass in the constant M/L ratio (de Vaucouleurs) com- 
ponent (solid curve with square points). The time delay values 
predicted by each lens model are plotted with triangles and con- 
nected by a dashed curve. These points correspond to the time 
delay scale on the right axis. A modest trend favoring dark-matter 
dominated models is visible. 

system ([Morgan et al.ll2008D in Pig. S While our new 
measurement is statistically consistent with our previ- 
ous estimate, the new expectation value is nearly 0.6 dex 
larger than our original analysis and has nearly twice 
the precision. We attribute this difference to two ef- 
fects. First, while we failed to measure the time delay, 
the new data substantially reduced the allowed range to 
—30 days < AtAB < days, leading to a significant nar- 
rowing of the source size probability density. Flux ra- 
tios estimated for an incorrect delay have extra power on 
short time scales which the model tries to fit using a more 
compact source size, leading to a broadening of the source 
size distribution toward smaller sizes. Since the time de- 
lay signal with our four new seasons of monitoring data 
has effectively ruled out positive delays, there are fewer 
cases in which the Monte Carlo code successfully repro- 
duced the noise from an incorrect delay using a small 
accretion disk size. Second, the simultaneous analysis 
of the (admittedly) sparse X-ray light curves provides a 
very strong constraint on the effective velocity distribu- 
tion v e , which, in turn, constrains r s . Blac kburne et aLI 
(2012) gained similar statistical leverage from the use of 
X-ray photometry in their analysis. 

We have made significant improvements in our un- 
derstanding of the QJ0158-4325 continuum source, but 
we are still unable to measure the time delay. This 
may be due to inadequacies in our data or techniques, 
such as uncertainties in our photometric measurements 
and/or insufficient signal from the time delay com- 
pared to the signal from microlensing, or it could pos- 
sibly be due to our oversimpl ification of the accre- 
tion d isk model. For example, lAbolmasov fc Shakural 
(2012) demonstrate that general relativistic effects be- 
come important during caustic crossing events and can 
explain high-magnification features in the light curves 



of SBS 1520+530 and Q 2237+0305. Several such fea- 
tures are present in the light curves of Q J0158-4325 
and are not perfe ctly fit by our present analysis . To be 
fair, the results of lAbolmasov fc Shakural ()2012D require 
disks with large inclination angles, and we see no evi- 
dence for this in Q J0158-4325. Nevertheless, we are 
left to conclude that the modeling of high-amplitude mi- 
crolensing events may be improved by the inclusion of 
a fully relativistic disk model or a more sophisticated 
model for the motions of sta rs in the lens galaxy (e.g. 
IPoindexter fc Kochanekll2010D . 

At X-ray wavelengths, a consistent picture of the struc- 
ture of the quasar continuum emission region is begin- 
ning to emerge. In the four systems with quantitatively 
measured continu u m em ission regions, HE 0435-1223 
(Blac kburne et al.l [20ll. RXJ 1131-1231 (|Dai et al.l 
120101) . PG 1115+080 (rMorgan et al.l [20081) and now 
QJ 0158-4325, the X-ray continuum emission source is 
located at ~ 10 r g . Improving on our previous work, 
we performed the second complete analysis of the hard 
and soft X-ray variability, leading to independent esti- 
mates for the sizes of the soft and hard X-ray sources. 
While our hard and soft size estimates are formally con- 
sistent with each other, there are some indications that 
the hard X-ray source is smaller th an the soft source. As 
we discussed in iChen et al.l (|2012| ). these data were not 
designed to produce good statistical estimates for the 
hard and soft band separately, so we expect that adding 
more X-ray epochs with better photon statistics should 
lead to markedly improved X-ray size estimates. 
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